Condensation of 'composite bosons' in a rotating BEC 



1—5 



N.K. Wilkin (i^^) and J.M.F. Gunn (i^^) 

'^^ School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom. 

''^'European Synchrotron Radiation Facility, BP. 220, 38043 Grenoble Cedex 9, France 

(February 1, 2008) 

We provide evidence for several novel phases in the dilute limit of rotating BECs. By exact cal- 
culation of wavefunctions and energies for small numbers of particles, we show that the states near 
integer angular momentum per particle are best considered condensates of composite entities, in- 
volving vortices and atoms. We are led to this result by explicit comparison with a description purely 
in terms of vortices. Several parallels with the fractional quantum Hall effect emerge, including the 
presence of the Pfafhan state. 
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In rotating superfluid ^He a vortex lattice forms which, 
on scales large compared to the vortex lattice parameter, 
has a velocity field indistinguishable from a rigid body 
co-rotating with the container. The vortices only per- 
turb the fluid density significantly over a region of order 
the coherence length, ^, around the core of each vortex 
(of the order of an Angstrom). Hence the arrangement 
of the vortex lattice is governed by minimising the ki- 
netic energy of the fluid in the rotating frame. One may 
say that the potential energy is 'quenched' by the incom- 
pressibility of the fluid. 

In the Bose condensed alkali gases although so far it 
has proved difficult experimentally to investigate the ro- 
tational properties of the condensates, there has been a 
vigorous theoretical debate |jl|-^ about the stability (or 
otherwise) of vortices in the condensates. At a mean 
field level (appropriate for moderate density), the inho- 
mogeneity of the condensate density and the existence of 
surface waves due to the harmonic well makes the descrip- 
tion difficult. Nevertheless the interparticle potential en- 
ergy is still largely unaffected by the presence of vortices 
in the limit where the coherence length is small compared 
to the extent of the condensate: it is the kinetic energy 
(and the single-particle trap potential) which determines 
the vortex positions. 

In this Letter we show that when the coherence length 
is comparable to the extent of the condensate, completely 
new effects occur. This is due to the kinetic (and single 
particle trap) energy being quenched, by a combination of 
spherical symmetry and the special properties of the har- 
monic well. Hence the ground state in the rotating frame 
is determined by the interparticle interactions alone, rem- 
iniscent of the fractional quantum Hall effect. Indeed we 
find stable states that are related to those found in the 
Hall effect (albeit in the less familiar regime of filling frac- 
tion, 1/ > 1). These include 'condensates' of composite 
bosons of the atoms attached to an integral number of 
quanta of angular momenta, as well as the Laughlin and 
Pfaffian [Q states. 

In a rotating reference frame, the standard Hamilto- 
nian for N weakly interacting atoms in a trap is [g : 
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where we have used the trap energy, fi^J K/m — hujQ as 
the unit of energy and the extent, {h^ /MKY^^, of the 
harmonic oscillator ground state as unit of length. (Here 
M is the mass of an atom and K the spring constant 
of the harmonic trap.) The coupling constant is defined 
as 7] = inna^h^ /M K)~^^'^ where n is the average atomic 
density and a the scattering length. The angular velocity 
of the trap, w, is measured in units of the trap frequency. 

In the dilute limit 77 ^ 1, which implies in existing ex- 
perimental traps that the number of atoms, iV, would be 
10 ^ N ^ 1000. Then the average coherence (or healing) 
length is ^ ~ Xjy/fM — > 00. It has been shown previ- 
ously Pi that in this limit, the problem becomes two- 
dimensional and the Hilbert space may be truncated to 
the 'lowest Landau level' states 0, V-'m(^) c«: z™e~l^l /^, 
where m > and z = x -I- iy in the plane normal to ui. 
Indeed at a; = 1 the problem is identical to the Quantum 
Hall problem, with a; replacing the magnetic field. 

We have determined the exact ground state, its en- 
ergy, Eq(ll})^ and excitation gap, A, for the Hamiltonian 
Eq. (|l|) using Mathematica for TV < 8 and w < 1. In ad- 
dition we have determined numerically the lowest eigen- 
values for iV < 10 as a function of ld. Lq{uj), the angular 
momentum of the ground state, is plotted in Fig. ^ for 
A^ = 6 with 77= 1/A. Angular momentum remains a good 
quantum number as we have made no symmetry breaking 
ansatz. 

A corresponding plot for *He in a rotating container 
would show jumps in the expectation value of Lo{ll)) as 
successive vortices enter the system. The inhomogeneous 
density of the condensate in a trap leads to more com- 
plex, but similar, behaviour in a mean field treatment |q| 
(appropriate in the high density limit). There are a num- 
ber of important features in Fig. n^, which are common to 
all values of N which we have studied. 

Firstly, at L — N there is a state which corresponds B 
to one vortex. 
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where Zc = (X]i=i ^i)/-^ i^ the centre of mass eoordinate 
'. From this point we will omit nor- 
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and 

malisation factors and the ubiquitous e^l^l /^. This state 
has an interparticle interaction energy E = riN{N—2)/A 
and becomes stable at wi = (1 — A^77/4). (In addition at 
uji all N > L>1 states are metastable M.) 
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FIG. 1. Stable states for TV = 6 in the rotating frame with 

ri=l/N 

Almost all of the other stable states can be labelled by 
L = n{N — m) where n and m are nonnegative integers 
(this includes the Laughlin state n = N, m = I). These 
values are close to 'n- vortex states' {L = nN), a possibil- 
ity we will return to. However the actual wavefunctions 
for these states most closely resemble some ||l^ used in 
the theory of composite fermions [O in the Hall effect. 

We define 
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(Note: ^P"^-^ = UtlQN-liz^) and i^^^ = ntiQi(^O-) 
Then the states of high overlap with the true states at 
L = n{N—m) may be written as: 

N 
1pn,m[\Zi\) = / ^ Qn\Zj-^)Qn\Zj^) ■ ■ ■ Qn\Zj^j^_^^) 

jl<J2<---<j(N-m) 

Table 1 shows the overlaps of ipn,m with the true ground 
states for those L. Their construction ensures that angu- 
lar momentum is used economically to lower the energy: 
any given particle pair, i and j, will at most be associated 
with two factors of (zi — Zj). 

The interpretation of the states at L — n{N—m) is that 
a particle in association with n quanta of angular mo- 
mentum is a particularly stable entity in the vicinity of 



L = nN . As L is reduced N—m particles remain with all n 
quanta and m have all the angular momentum removed. 
This is our main result, which occurs at small accessi- 
ble angular momenta. This is reminiscent of the 'bound 
state' composite fermions of electrons and vortices [O. 
We will return to this point. 

We will now attempt to reconcile the composite bo- 
son states to the vortex states found in the Nonlinear 
Schrodinger equation 0]. The following argument indi- 
cates a connection. Consider incompressible irrotational 
fluid ('Helium') in a two-dimensional circular container, 
of radius i?, with n point vortices at radial coordinates 
Tq, . There the angular momentum of the fluid is [nsl : 
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Le. the angular momentum is reduced from L = nN by 
the vortices being off-centre. 

To test this notion, we firstly localise the vortices (re- 
sulting in a non-rotationally invariant state) by superpos- 
ing states with different L. Using L — IQ and L = 8 for 
N — 5 (the Pfafhan state rules out A^ = 6) yields the con- 
tour plot of probability density, Fig.0. The two dimples 
might be interpreted as two off-centre vortices (hence the 
angular momentum is lower than L = 2N). The figure is 
reminiscent of the figures in |8J , although the changes in 
density are rather small by comparison. Note however, 
the superposition is certain to create features periodic 
with cos 26, where 9 is the polar angle. 
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FIG. 2. The probability density for the superposition of 
the states L = 2N and L = 2N -2 for TV = 5 with 'dimples' 
reminiscent of vortex cores. 



To quantify these ideas, we introduce the lowest Lan- 
dau level vortex factors, 0^=1 (Ca ~ ^i) with complex vor- 
tex coordinate C,^- A complete set of particle states of 
angular momentum L is obtained using n = L vortices: 
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where \C,\^ = Xla=i ICaP- This follows from the result 
p4| that an nth order polynomial in the z's may be ex- 
panded using the products of powers of the elementary 



symmetric polynomials, Cr, < r < n: 
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Although there is no need to write wavefunctions of both 
the particle, Zi, and the vortex coordinates, Cq, it will be 
convenient. 

It will be useful to note the form of the resulting par- 
ticle wave function, ip, when there is one vortex. The 
construction implies the natural one-vortex states are 
'f'piC*) = C*^/(p!^)- The corresponding particle states 
are Viv-p fx Cp{{zi}) ior < p < N; and for p > N 
ipN-p = 0. Note that the angular momentum of the par- 
ticles is L = N — p, consistent with the special cases: 
TpN = 1 i.e. placing a vortex in this state has no effect on 

the non-rotating condensate; ^n-i =^ Si=i ^i ~ ■^c^ ^^'^ 
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i.e. a 'simple' single vortex. 



tcx factors (even more for larger L). In addition, the 
number of vortices is not fixed, as the number in the 
vortex state (^ is indeterminate since they do not af- 
fect the particle wave-function (in a sense it is the vortex 
vacuum state). This is in stark contrast to the incom- 
pressible (^ -^ 0) case where the number of vortices is 
fixed and they are classical entities. 

It might be supposed that although there may be a 
fluctuating vortex population at large vortex quantum 
numbers, this is in the tail of the particle wavefunction 
and the description may be simple near the centre of the 
trap. This is determined by computing the single- wortea; 
density matrix, p^{C,C); for the exact state L = 2iV— 2, 
7V = 6, Fig. I 
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FIG. 3. Eigenvalues (pm) of the single vortex density ma- 
trix for eigenfunctions C" for L = 27V— 2 for N = 6. The trace is 
normalised to unity, having suppressed the weight associated 
with the vortex 'vacuum' state, m = N. 

If the displaced vortex picture were correct, one might 
expect a factor in the vortex wavefunction of the form 



It is tempting to relate the stable states at L = 
n{N — m), to Eq. 0, as being n vortices in state p = m, 
with associated displacement of the vortex determined 
by (ICP)p = P + 1. However this implies r^ — p + 1 and 
this leads to a contradiction unless L ^ N'^ (using the 
empirical relation p = m < n), which is too restrictive. 

Moreover, this purely vortex description, 0(C), requires 
more vortices than n [in L = n{N~m)). For example 
L = N: expanding the product in Eq|| we see there is a 
term z^ = Ci{{zi})^ whose generation requires N vor- (Noting e 



(Ci ~ Q)'^ corresponding to the vortices rotating around 
the centre of the trap. This factor alone would lead to 
the following eigenvalues, p^ (with corresponding eigen- 
vectors C"), of p^: pl^= I, pl = ^ and p^ = \. As 
can be seen from Fig.ra, this is not the case. The most 
pronounced feature is a maximum at ?7i = 0. The vortices 
tend to condense, in the m — Q state, not to separate in 
ICr ~ Q.\- (Further evidence comes from evaluating the 
particle density matrix [|j). 

These difficulties in describing the V'n,m states purely 
in vortex variables occur because the particles are binding 
to the vortices. This leads to a strongly correlated state 
whose description requires additional vortex variables if 
they are used alone. One interpretation uses ideas from 
the quantum Hall effect [Q : At the centre of each vortex 
there is a decrease in the particle density. Thus in terms 
of interparticle interactions, this is a low energy region 
for an additional particle. 

Mathematically this is described most easily for the 
Laughlin state, using N vortices, Coi with a factor 
Jljj /Q,(Ca — Zi) where the a-th particle experiences no 
suppression of its amplitude: it is 'bound'. This can be 
expressed as 
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e^^ respectively play the roles 
]l5|,p| of a delta function and its nth derivative within the 
lowest Landau level.) I.e. the delta function factors bind 
the ith particle to the ith vortex. 

To generate the states ijjn,m we use the derivatives of 
the Lowest Landau level delta function so that: 
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which can be interpreted as a 'condensate' of {N — m) 
composites (each consisting of an atom and a vortex) 
in the state, C,*"^ with r — [N—l — n). The remaining 
unbound atoms remain condensed in the single particle 
ground state. (The states L = n{N — m) are also selected 
using a composite fermion approach pq.) 

The remaining stable states are consistent with the 
bosonic Pfaffian state ^^, at L = ^N{N-2) for even 
TV and L = ^(^-1)^ for odd N. 

^'''i{z.})^Y[iZr~Z,)Pi 
i<j 

where the Pfaffian is defined 
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where A denotes antisymnietrisation of the following 
product. (This is generalised for odd N by omitting one 
of the particles in each term of the antisymmetrisation 
||14|.) The overlaps of ip^^ with the exact ground state 
for TV = 5 and L = 8, TV = 6 and L = 12 and A^ = 7 and 
L = 18 are: OJl^, OJO^ and O.BO^. 

Some nearby stable states, e.g. L — 10 and L = 14, 
are well described as simple modifications of the Pfaf- 
fian state. This uses the conjecture (which has been 
demonstrated by direct evaluation for 4 < TV < 8) that 
the PfafRan state may be represented by a product of 
two Laughlin states for N/2 particles (or for odd N, a 
cluster of (iV-l)/2 and one of (TV + l)/2): 
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where the two subsets, cti and tT2, each have N /2 particles 
((TV-l)/2 and (7V+l)/2 for odd TV). S indicates that the 
wave function is symmetrised over the distribution of the 
particles into these subsets. These two well-correlated 
clusters appear to be 'dual' to the clusters of Halperin 
p8[ which have a high internal energy, due to the lack of 
nodal factors. 

For example the state TV = 6, L = 14 has overlap 
0.96^ with a state with two quanta of angular momenta 
in the centre of mass motion of the clusters (defining 
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1 or 2): 



2„;,Pf 



^'^'%{z^}) - s{z, - z^Yi^ 



The state TV = 6, L = 10 has overlap 0.97^ with a state 
where there is one factor of centre of mass motion and 
one vortex has been 'removed' from one of the clusters: 
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(The apparent asymmetry of the last factor only involv- 
ing the first cluster, cti, is illusory due to the overall sym- 
metrisation.) 

In conclusion, this Letter provides evidence that the 
weak coupling limit of rotating EEC's contains some 
novel phenomena. These occur even in the regime where 
one would anticipate small numbers of vortices and hence 
should be open to experimental investigation in the near 
future. 



ACKNOWLEDGMENTS 



We would like to thank J.T. Chalker, N.R. Cooper, 
D.S. Rokhsar and R. A. Smith for helpful discussions. 
We are grateful to P. Carra and the ESRF for hospitality 
whilst this work was completed and ITP (PHY-94-07194) 
during the early stages of the work. We acknowledge fi- 
nancial support from the EPSRC GR/L28784,GR/K6835 
and GR/L29156. 





2N-i 


27V-2 


2N 


3iV-6 


3Af-3 


3N 


47V-8 


iN-i 


AN 


5W-10 


5Af-5 


5N 


QN-e 


GN 


5 




8 Ilk 

0.98^ 


10 

0.872 




12 

0.992 


15 

0.962 






20 

1 












6 




10 4k 

0.86^ 


12 4k 

0.692 


12 4k 

0.792 


15 

0.952 






20 

0.992 


24 

0.942 






30 

1 






7 




12 

0.832 




15 

0.492 


18 4k 

0.852 






24 






30 


35 




42 

1 


8 


12 

0.662 


14 




18 




24 ♦ 


24 4k 




32 


30 


35 




42 





TABLE 1. Stable states for TV < 8: the upper number is their angular momentum and the lower is their overlap with the Q 
wavefunctions. tfi indicates that the wavefunction can also be written (or derived from) a Pfaffian state. 
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